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The relaxation of a dewetting contact line is investigated theoretically in the so-called 
"Landau-Levich" geometry in which a vertical solid plate is withdrawn from a bath of 
partially wetting liquid. The study is performed in the framework of lubrication theory, in 
which the hydrodynamics is resolved at all length scales (from molecular to macroscopic) . 
We investigate the bifurcation diagram for unperturbed contact lines, which turns out 
to be more complex than expected from simplified 'quasi-static' theories based upon an 
apparent contact angle. Linear stability analysis reveals that below the critical capillary 
number of entrainment, Cac, the contact line is linearly stable at all wavenumbers. Away 
from the critical point the dispersion relation has an asymptotic behaviour a oc \q\ 
and compares well to a quasi-static approach. Approaching Cac, however, a different 
mechanism takes over and the dispersion evolves from ^ |g| to the more common . 
These findings imply that contact lines can not be treated as universal objects governed 
by some effective law for the macroscopic contact angle, but viscous effects have to be 
treated explicitly. 



1. Introduction 

Wetting and dewetting phenomena are encountered in a variety of environmental and 
technological contexts, ranging from the treatment of plants to oil-recovery and coating. 
Yet, their dynamics can not be captured within the framework of classical hydrodynamics 
- with the usual no-slip boundary condition on the substrate - since the viscous stress 
diverges at the contact line (|Huh fc Scriven 19711 IDussan et al. 1974^ . The description 
of moving contact lines has remained a great challenge, especially because it involves 
a wide range of length scales. In between molecular and millimetric scales, the strong 
viscous stresses are balanced by capillary forces. In this zone, the slope of the free 
surface varies logarithmically with the distance to the contact line so that the interface 
is strongly curved, even down to small scales (jVoinov 1976llCox 1986p . Ultimately, the 
intermolecular forces due to the substrate introduce the physical mechanism that cuts 
off this singular tendenc y (IVoinov 19761 ICox 19861 Ide Gennes 19 86 Blak e et al. 1995|. 
IPismen fc Pomeau 2000|l . 

A popular theoretical approach has been to assume that all viscous dissipation is local- 
ized at the contact line, so that macroscopically the problem reduces to that of a static 
interface that minimizes the free energy. In such a quasi-static approximation one does 
not have to deal explicitly with the contact line singularity: the dynamics is entirely gov- 
erned by an apparent contact angle 9a that serves as a boundary condition for the inter- 
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Figure 1. (a) A standard geometry to study contact line dynamics is that of a vertical solid 
plate withdrawn from a bath of liquid with a constant velocity U. The position of the contact 
line is indicated by z^i- (b) In this paper we study the relaxation of transverse perturbations of 
contact lines, by computing the evolution of the interface profile h{z,y,t). 



face (jVoinov 19761 ICox 19861 |Joanny fc de Gennes T9841 [Goiestanian fc Raphael 2001bl 
[Nikolayev fc Beysens 2003 1. This angle is a function of the capillary number Ca = r]U/j, 
which compares the contact line velocity U to the capillary velocity "f/rj, where 7 and 
77 denote surface tension and viscosity. Since the dissipative stresses are assumed to 
be localized at the contact line, viscous effects will modify the force balance determin- 
ing the contact angle. This induces a shift with respect to the equilibrium value Oe- 
Within this approximation, the difficultly of the contact line problem is hidden in the 
relation 9a{Ca), which depends on the mechanism releasing the singularity. While it 
is agreed upon that the angle increases with Ca for advancing contact lines and de- 
creases in the receding case, there are many different theories for the explicit form 
(IVoinov 19761ICOX 19861 Ide Gennes 19861 IBlake et al. 19951) . 

Experimentally, however, it has turned out to be very difficult to discriminate between 
the various theoretical proposals (Hoff maii 19751 ILe Grand et al. 2 005 Ri o et al. 2005p . 
All models predict a nearly linear scaling of the contact angle in a large range of Ca and 
the prefactor is effectively an adjustable parameter (namely the logarithm of the ratio 
between a molecular and macroscopic length). Differences become more pronounced 
close to the so-called forced wetting transition: it is well known that the motion of 
receding contact lines is limited by a maximum speed beyond which liquid deposition 
occurs (jBlake fc Ruschak 1979[ Ide Gennes 1986|l . An example of this effect is provided 
by drops shding down a window. At high velocities, these develop singular cusp-hke tails 
that can emit little droplets ( [Podgorski et al. 200H ILe Grand et al. 2005p . Similarly, 
solid objects can be coated by a non-wetting liquid when withdrawn fast enough from 
a liquid bath (jBlake fc Ruschak 19791 |Quere 19 91 Scdc v fc Petrov 199ip , see Fig. [l}i. 
Above the transition, a capillary ridge develops (Snoeijer et al. 20061) that eventually 
leaves a Landau-Levich film of uniform thickness (jLandau & Lcvich 1942|. 

An important question is to what extent a quasi-static approximation, in which dissi- 
pative effects are taken localized at the contact line, are able to describe these phenom- 
ena. Only recently, the problem has been addressed by using a fully hydrodynamic model 
that properly incorporates viscous effects at all length scales ( [Hocking 2001|[Eggers 2004} 
[Eggers 2005| ). It was found that stationary meniscus solutions cease to exist above a crit- 
ical value Cac, due to a matching problem at both ends of the scale range: the highly 
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Figure 2. Experimental realization of contact line perturbations. The contact line is deformed 
by "wetting defects" on the partially wetting plate. The narrow connection between the de- 
fect and the bath undergoes a Rayleigh-Plateau-like instability, leaving a periodically deformed 
contact line. 



curved contact line zone and the macroscopic flow fEgg ers 2004[ |Eggers 2005p . An in- 
teresting result of this work is that both the value of C&c and the emerging 9a (Ca) are 
not universal: these depend on the inclination at which the plate is withdrawn from the 
liquid reservoir. Hence, the large scale geometry of the interface does play a role and the 
dynamics of contact lines can not be captured by a single universal law for 6*0 (Ca). 

Golestanian and Raphael identified another sensitive test to discriminate contact line 
models ( Golestanian fc Raphael 2001a Golestanian fc Rap hael 2001b|). They considered 



the relaxation of dewetting contact lines perturbed at a well-defined wave number q, as 
shown in Fig. This can be achieved experimentally by introducing wetting defects 
on the solid plate, separated by a wavelength A ( [Ondarguhu fc Veyssie 199T| . As can be 
seen in Fig. |2] these defects create a nonlinear perturbation when passing through the 
contact line, but eventually the relaxation occurs along the Fourier mode with q = 2tt/X 
()Delon et al. 2006^ . Using a quasi-static theory, Golestanian and Raphael predict that 
the perturbations decay exponentially ^ e~'^* with a relaxation rate 

a=|g|^/(Ca), (1.1) 
V 

where /(Ca) is very sensitive to the form of Oa{Ga). Their theory is built upon the 
work by |Joanny fc de Gennes 1984) who already identified the scaling proportional to 
\q\ for static contact lines (Ca = 0). Ondarguhu and Veyssie experimentally confirmed 
this Igl dependence in the limit of Ca — QOndarguhu fc Veyssie 199T| , while more 
recently, it has been argued that this scaling should saturate to the inverse capillary 
length ly — y^^/pg in the large wavelength limit ( [Nikolayev fc Beysens 2003| . However, 
an intriguing and untested prediction for the dynamic problem is that the relaxation 
times diverge when approaching forced wetting transition at Cac. This "critical" behavior 
should occur at all length scales and is encountered in the prefactor /(Ca), which vanishes 
as Ca ^ Ca^. 

In this paper we perform a fully hydrodynamic analysis of perturbed menisci when a 
vertical plate is withdrawn from a bath of hquid with a velocity U (Fig. [T]) . Using the 
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lubrication approximation it is possible to take into account the viscous dissipation at 
all length scales, from molecular (i.e. the slip length) to macroscopic. We thus drop the 
assumptions of quasi-static theories and describe the full hydrodynamics of the problem. 
The first step is to compute the unperturbed meniscus profiles as a function of the plate 
velocity. We show that these basic solutions undergo a remarkable series of bifurcations 
that link the effect that stationary cease to exist beyond Cae dEggers 2004| ), to the recently 
observed upward propagating fronts beyond the transition (Snoeijer et al. 2006). Then 
we study the dispersion of contact line perturbations through a linear stability analysis. 
Our main findings are: (i) the relaxation time for the mode q = scales as |Ca— Cad"^/^; 
(ii) finite wavelength perturbations always decay in a finite time even right at the critical 
point; (iii) the scaling a cx |g| proposed by Eq. (1.1 1 breaks down when approaching C&c- 
These results illustrate the limitations of simplified theories based upon an apparent 
contact angle. 

The paper is organized as follows. In Sec. [2] we summarize the results from a quasi- 
static theory and generalize the work by Golestanian & Raphael. The heart of the 
paper starts in Sec. |3] where we formulate the hydrodynamic approach and compute the 
bifurcation diagram of the base solutions. After addressing technical points of the linear 
stability analysis in Sec. |4] we present our numerical results for the dispersion relation 
in Sec. [5j The paper closes with a discussion in Sec. [6] 



2. Results from quasi-static theory 

We briefiy revisit the quasi-static approach to contact line perturbations, which will 
serve as a benchmark for the full hydrodynamic calculation starting in Sec.|3] The results 
below are based upon the analysis of Golestanian & Raphael, which has been extended 
to long wavelengths and large contact angles. 

2.1. Short wavelengths: ql-y S> 1 

At distances well below the capillary length, Ij — ^/j/pg, we can treat the unperturbed 
profiles as a straight wedge of angle 0a- Perturbations should not affect the total Laplace 
pressure, and hence not the total curvature of the free interface. The interface will thus 
be deformed as sketched in Fig. [3^: the advanced part of the contact line has a smaller 
apparent contact angle than the unperturbed 0a- According to ^^(Ca), such a smaller 
angle corresponds to a higher velocity with respect to the plate, hence the perturbation 
will decay. From this argument one readily understands that the rate of relaxation a, 
depends on how a variation of induces a variation of Ca, and thus involves the derivative 
dCa./d0a dColestanian Ik Raphael 2003p . 

Working out the mathematics, see Appendix |B] we find 

■qa _ tangg / dtan^g X^^ 

1^ ^ ^cos0a V ^^Ca ) ' ^'^ 
This implies that the timescale for the relaxation is set by the length q^^ and the capillary 
velocity "f/rj, where 7 represents surface tension and rj is the viscosity. We therefore 
introduce 

aoo(Ca)= hm , (2.2) 

which will be used later on to compare to the hydrodynamic calculation in the limit of 
large g. 
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Figure 3. (a) Macroscopic representation of the interface shape near a perturbed contact hne. 
The advanced part of the contact hne has a smaller apparent contact angle than the unperturbed 
6a, and will thus have a higher speed with respect to the plate. This will decrease the amplitude 
of the perturbation, (b) Cross-sections of the perturbed interface profile along z. Note that the 
interface joins the static bath at z = 0. 



2.2. Large wavelengths: qlj 1 

When considering modulations of the contact line with 1 /q of the order of the capillary 
length, one can no longer treat the basic profile as a simple wedge. Instead, one has to 
invoke the full profile /io(-2) and the results for a are no longer geometry independent (see 
also ISekimoto et al. 1987P . For the geometry of a vertical plate immersed in a bath of 
liquid, we can characterize the profiles by the "meniscus rise", indicating the position of 
the contact line Zd above the liquid bath (Fig.jlJ. This is directly related to the contact 
angle as (jLandau fc: Lifschitz 1959^ 



^ ±1^^2(1 -sm9a), (2.3) 

where the sign depends on whether 9^ < tt/2 (positive), 9a > 7r/2 (negative). In fact, 
this relation is often used to experimentally determine ^^(Ca), since the meniscus rise 
Zc\{Ca,) can be measured more easily than the slope of the interface. 

We now consider the relaxation rate cto for perturbations with q = 0. Such a per- 
turbation corresponds to a uniform translation of the contact line with Az. Using the 
empirical relation Zci(Ca), we can directly write 

-1 

Az. (2.4) 



dAz 



1 



dt 



= -- ACa 
V 



--if .^^Y 

r] \dCa.J 



Hence, 

7 

In terms of the contact angle, using Eq. 



I2T3I), this becomes 



sin 9a) f dtan 9, 



(2.5) 



(2.6) 



7 I cos-^ 9a I V '^Ca 

Besides some geometric factors, this result has the same structure as the relaxation for 



small wavelengths, Eq. (2.11. The crucial difference, however, is that the length scale of 



the problem is now L instead of q. Comparing the relaxation of finite wavelenghts, tT„ 
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with the zero mode relaxation, ctqi we thus find 



(J. 



ql^g{Oa) for > 1, (2.7) 



where the prefactor g{Oa) reads 



= i ;^"^^"!^'"^- . (2.8) 

V2(l~sin0J 



Using the definition Eq. (2.2) we thus find the quasi-static prediction 



= aogi9a). (2.9) 



2.3. Physical implications and predictions 



The predictions of the quasi-static approach can be summarized by Eqs. (2.11, (2.51, 



and (2.9 1. The relation cr cx \q\ was already found by [Joanny fc de Gennes 1984[ who 
referred to this as the "anomalous elasticity" of contact lines. The linear dependence on 
q contrasts with the more generic scaling q^ for diffusive systems, and has been confirmed 



experimentally by |Ondarguhu fc Veyssie 1991 in the static limit, Ca = 0. An interesting 



consequence is that the Green's function corresponding to this dispersion relation is a 
Lorentzian cx (1 -t- [y /w{t)]'^)~^ , whose width w{t) grows linearly in time. The prediction 
is thus that a localized deformation of the contact line, similar to Fig. [2] but now for 
a single defect, will display a broad power- law decay along y. In the hydrodynamic 
calculation below we will identify a breakdown of this phenomenology in the vicinity of 
the critical point. 

On the level of the speed-angle law ^^(Ca), the wetting transition manifests itself 



through a maximum possible value of Ca, i.e. dO/dCa, = oo. According to Eqs. (2.11 



and (2.6 1, this suggests a diverging relaxation time ct ^ at all length scales. Assume that 
the scaling close to the maximum is 9a — 9c (Cac — Ca)'', where generically one would 



expect P = 1/2. If the critical point occurs at zero contact angle, 9c = 0, Eq. (2.1 1 yields 
a scaling dq cx (Cac — Ca) for the case of large q. For q — or when 9c ^ 0, one finds 
Oq cx (Cac — Ca)^"''. Below we show that the mode g = indeed displays the latter 
scaling with /3 = 1/2. However, the relaxation times of finite q perturbations always 
remain finite according to the full hydrodynamic calculation, even at the critical point. 



3. Hydrodynamic theory: the basic profile /lo(^) 

This section describes the hydrodynamic theory that is used to study the relaxation 
problem. After presenting the governing equations, we reveal the nontrivial bifurcation 
diagram of the stationary solutions, h^i^z), for different values of the capillary number. 
This allows an explicit connection between the work on the existence of stationary menisci 
( [Hocking 200l][Eggers 2004l|Eggers 2005 ), and recently observed transient states in the 



deposition of the Landau-Levich film ( Snoeijer et al. 2006] ). All results presented below 



have been obtained through numerical resolution of the hydrodynamic equations using 
a Runge-Kutta integration method. 

3.1. The lubrication approximation 

We consider the coordinate system (z, y) attached to the solid plate, as indicated in 
Fig.[l] The position of the liquid/vapor interface is denoted by the distance from the plate 
h{z,y,t). To really cover the range of length scales from molecular to millimetric, the 
standard approach is to describe the hydrodynamics using the lubrication approximation 
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(lOron et al. 1997|) . This is a long wavelength expansion of the Stokes flow based upon 
Ca <C 1, which reduces the free boundary problem to a single partial differential equation 
for h{z,y,t). Of course, one still has to deal explicitly with the fact that viscous forces 
tend to diverge as h —> 0. Here we resolve the singularity by introducing a Navier slip 
boundary condition at the plate, 

Vz = Is ^ , (3.1) 

that is characterized by a slip length Ig . Such a slip law has been confirmed experimen- 
tally, yielding values for Is ranging from a single molecular length up to a micron depend- 
ing on wetting properties of the liquid and roughness of the solid (Barrat & Bocquet 1999 



ICottin-Bizonne et al. 2005llPit et al. 2000,, Thompson fc Robbins 1989). Different mech- 



anisms releasing the contact line singularity will lead to similar qualitative results, as long 
as the microscopic and macroscopic lengths remain well separated. 

The lubrication equation with slip boundary condition reads (jOron et al. 1997]l 

dth + W ■ {hlJ) = , (3.2) 
_ 3r/({7e^-U) 

Here U is the plate velocity, U(z, y, t) = UzGz + UyGy is the depth-averaged fluid velocity 
inside the film, while V = ezdz + Sydy. The first equation is mass conservation, while the 
second represents the force balance between surface tension 7, gravity pg, and viscosity 
r], respectively. We maintain the full curvature expression expression 

_ {l + dyh"^) dzzh+ {1 + dzh'^) dyyh-2dyhdzhdyzh 

" {1 + dzh^ + dyh^)y^ ' ^ ' 

which allows a proper matching to the liquid reservoir away from the contact line. 

In the remainder we rescale all lengths by the capillary length 1^ — y''j/ pg, and all 
velocities by j/rj, yielding the dimensionless equations 

dth + W ■{hV) = 0, (3.5) 

, 3(Cae,-U) 
VK-e^H — ^, , = 0. (3.6) 

h{h + 3ls) 

The timescale in this equation thus becomes rjl^/^. 

3.2. Boundary conditions 

We now have to specify boundary conditions at the liquid reservoir and at the contact 
line - see Fig. [3j Far away from the plate, h ^ 00, the free surface of the bath is 
unperturbed by the contact line. Defining the vertical position of the bath at z = 0, we 
can thus impose the asymptotic boundary conditions as z — s- 

dzh = —00, 
dyh = 0, 

K = 0. (3.7) 

We impose at the contact line, at z — z^i, that 

\Vh\ = tan^ci, 

hV = 0. (3.8) 
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The first condition determines the position of the contact line, while the third condition 
ensures that no liquid passes the contact line. This condition is not trivial, since the equa- 
tions admit solutions where the liquid velocity diverges as ~ l/h. The second condition 
imposes the microscopic contact angle, 9ci, emerging from the force balance at the contact 
line. This condition is actually hotly debated: for simplicity it is often assumed that this 
microscopic angle remains fixed at its equilibrium value ( [Hocking 2001 Eggers 2004^ , 



but measurements have suggested that this angle varies with Ca (jRame et al. 2004| . We 
will limit ourselves to presenting an argument in favour of a fixed microscopic angle. The 
boundary condition arises at a molecular scale, ZvdWi which the fiuid starts to feel 
the van der Waals forces exerted by the substrate. This effect can be incorporated by a 
disjoining pressure A/h^, where the Hamaker constant A (x '^ly^-^y (IsnMslachvili 1992^. 



At h — lydw, this yields a contribution of the order Ah'/ly^T^ in Eq. (3.3). Taking 
^Vdw ~ hi the viscous stresses will have a relative infiuence on the disjoining term, and 
thus on the boundary condition, of the order of Ca ^ 10"^, so that the contact angle 
should roughly remain within 1% of its equilibrium value. This analysis of the micro- 
scopic contact angle will be extended and compared to novel experimental results in a 
forthcoming paper Pelon et al. 2006|l . 

3.3. The basic profile and the bifurcation diagram 

We first solve for the basic profile hQ{z), corresponding to a stationary meniscus that 
is invariant along y. From continuity we find that hUz is constant, which using the 



boundary condition Eq. (3.8 1 yields Uz = 0. The momentum balance Eq. (3.6 1 for HqIz) 
thus reduces to 

with 

^^0= (3.10) 

Close to the bath, i.e. z ~ 0, the height of the interface becomes much larger than the 
capillary length, so that one can ignore the viscous term. The asymptotic solution of 



Eq. (3.9 1 near z — thus simply corresponds to that of a static bath. 



/iQ = — In z/c, 
1 



Kq = Z, (3-11) 



which indeed respects the boundary conditions of Eq. (3.7 1. This solution has one free 
parameter, c, that can be adjusted to fuUfiU the remaining boundary condition at the 
contact line, Hq — — tan^ci- It is not entirely trivial that this uniquely fixes the solution 



since Eq. (3.9 1 degenerates for Hq = 0. One can show that the asymptotic solution for 



small Z = Zc\ — z reads ( [Buckingham et al. 2003[ ) 

hQ^t^ne,,Z-^^^^^^^^^4l^^ZHnZ+~cZ\ (3.12) 
2ls tan ^ci 

where c is the sole degree of freedom. The expansion can be continued to arbitrary order 
once Oci and c are fixed. 

So, the basic profile ho{z) is indeed determined by the microscopic parameters 9ci and 
Is, and the (experimental) control parameter Ca. This is illustrated i n Fig. [4^, show ing 
Zci(Ca) for fixed parameters ^ci — 51.5° and h — 5 ■ 10"''. Similar to [Hocking 2001 [ and 



Relaxation of a dewetting contact line Part 1: A full-scale hydrodynamic calculation^ 
(a) 5^ (b) I 




T \ \ 1 

210-3 410-3 610-3^ 810-3 

Ca 




— I 1 1 1 

210-3 410-3 610-3^ 810-3 

Ca 



Figure 4. (a) Contact line position z^i at equilibrium as a function of the capillary number Ca 
for fixed parameters 6c\ = 51.5° and Zs = 5 • 10~^. Stationary solutions cease to exist above a 
critical value Cac — 0.00759- • •. The horizontal bar denots z^i — \/2. (b) Same as (a), but now 
showing the full range of Zd- The solutions undergo a sequence of saddle-node bifurcations with 
ultimately Zd oo, with a corresponding Ca* = 0.00693- • • (dashed line). 



[Eggers 2004[ we find that stationary meniscus solutions only exist up to a critical value 
Cac- Beyond this capillary number the interface has to evolve dynamically and a Uquid 
film will be deposited onto the plate. One can use Fig.|4^ to extract the apparent contact 
angle 9a, via Eq. (2.3 1. The critical capillary number is attained when Zc\ = 1.4076- • •, 
which is very close to \/2 = 1.4142 • • •. This confirms the predictions by [Eggers 2004| 
that stationary solutions cease to exist at a zero apparent contact angle. This slight 
difference from is due to the fact that Eggers's asymptotic theory becomes exact 
only in the limit where Is — > 0, so that minor deviations can indeed be expected. Other 
values of 9ci and Is lead to very similar curves, always with a transition at Zd — V^, 
but with shifted values of Cac. This critical value roughly scales as Cac oc 9^i/lnlj^ 
dEggers 2004| [Eggers 2005| . 

In fact, the existence of a maximum capillary number is due to a saddle-node bifurca- 
tion, which originates from the coincidence of a stable and and unstable branch (this will 
be shown in more detail in Sec. [S]). As can be seen from Fig. |4|d, there is a branch that 
continues above Zd = V^- Surprisingly, these solutions subsequently undergo a series 
of saddle-node bifurcations, with capillary numbers oscillating around a new Ca*. This 
asymptotically approaches a solution of an infinitely long flat film behind the contact 
line. Figure |5] shows the corresponding profiles ho, and illustrates the formation of the 
film. This film is very different from the so-called Landau-Levich film, which was com- 
puted in a classic paper ILandau fc Levich 1942|) . The Landau-Levich solution is much 
simpler in the sense that it does not involve a contact line and does not display the non- 
monotonic shape shown in Fig. [5] The difference markedly shows up in the thickness of 
the film: while the Landau-Levich film thickness scales as Ca^^^ , the film with a contact 
line has a thickness hoo = VsCa*. Note that very similar film solutions were already 
identified by [Hocking 2001 and more recently by IMiinch fc Evans 20051 in the context of 
Marangoni-driven flows. 

These new film solutions have indeed been recently experimentally, as transient states 
in the deposition of the Landau-Levich film ( |Snoeijer et al. 2006 p . In fact, the transition 
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Figure 5. Various basic solutions ^o(^) along the bifurcation diagram of Fig. H see inset. 



towards entrainment was observed to coincide at Ca* , hence well before the critical point 
Ca.c and with 6a ^ 0. For fibres, on the other hand, the condition of vanishing contact 
angle has been observed experimentally by ISedev fc Petrov 19911 We come back to this 
issue at the end of the paper. 



3.4. Physical meaning of the apparent contact angle 



From the definition of Eq. (2.3 1, it is clear that the apparent contact angle represents 
an extrapolation of the large scale profile using the static bath solution. In reality, 
however, the interface profile is strongly curved near the contact line and the contact 
angle increases to a much larger 9c\. This has been shown in Fig. [6] revealing the 
logarithmic evolution of the interface slope close to the contact line. This is different from 
the static bath solutions, for which the slope decreases monotonically when approaching 
the contact line (Fig. |6|d, dashed curve). Another way to define a typical contact angle in 
the dynamic situation could thus be to use the inflection point, which yields the minimum 
slope of the interface. However, when using 9a(Ca,) as an asymptotic matching condition 
for an outer scale solution, like in a quasi-static theory, it is clear that it only makes 
sense to use the extrapolated version. 



4. Linear stability within the hydrodynamical model 

We now turn to the actual linear stability analysis within the hydrodynamic model. 
This section poses the mathematical problem and addresses some technical issues related 
to asymptotic boundary conditions. The numerical results will be presented in Sec. |5] 

4.1. Linearized equation and boundary conditions 
We linearize Eqs. ( 3.5|3.6 1 about the basic profile /io(z), writing 

h{z, y, t) = /lo(z) + e h^{z) e--*+»92', (4.1) 

n{z, y, t) ^ Ko{z) + e k^{z) e-*+-iy ^ (4.2) 

f/,(z,y,i) (z)e-^*+^«^ (4.3) 

C/j,(z,y,i)=eC/j,i(z)e-"*+^9^ (4.4) 
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Figure 6. (a) Variation of the slope h'g as a function of the distance to the contact hne, for 
different values of Ca. (b) The profile for the critical solution (solid line), compared to the 
static bath with da ~ (dashed line). 



Here we used that the basic velocity Uq = 0, so that the velocity is of order e only. 
The linearized equation is not homogeneous in z, due to hg^z), so the eigenmodes are 
nontrivial in the z direction. From the y-component of Eq. (3.6 1, one can eliminate Uyi 
in terms of ki, as 



Uyi{z) = i iqho {ho + 3ls) Ki{z). 
It is convenient to introduce the variable 

Fiiz) = hoiz)U,i{z), 



(4.5) 



(4.6) 



which represents the flux in the z direction at order e (the zeroth order flux being zero) . 
Writing the vector 



X 



(4.7) 



Kl 

V ^1 J 

one can cast the linearized equation for the eigenmode as 

d,X = AX, (4.8) 

where dz denotes the derivative with respect to z and the right hand side is a simple 
matrix product. From linearization of Eqs. (3.4 1, (3.5 1, and (3.6 1 one finds 

/ 1 

3Ca i-Q 3ls \ r\ n 

\ a 



q^hliho + 3ls)/3 
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The eigenmodes and corresponding eigenvalues Uq of this hnear system are determined 
through the boundary conditions. At the contact hne, we have to obey the boundary 



conditions of a microscopic contact angle tan6'ci and a zero flux (see Eq. (3.8 1). To 
translate this in terms of the linearized variables, we have to evaluate | V/i| at the position 
of the contact line Zd+Az. Along the lines of App.[B]one flnds Az = ee~'^*+*'^ft.i/ tan 9c\. 
Linearizing |V/i| then yields the boundary conditions for the eigenmode 

^1 =^ tangei ^'-^^ 

Fi ^ 0. (4.10) 

At the side of the bath, z 0, the conditions become 

iqhi = Q q = V hi^O, (4.11) 

Kl = 0. (4.12) 

Below we identify the two relevant asymptotic behaviors at the bath respecting these 
boundary conditions. 



4.2. Shooting: asymptotic behaviors at bath 
The strategy of the numerical algorithm is to perform a shooting procedure from the bath 



to the contact line, where we have to obey the two conditions Eq. (4.9) and (4.10). We 
thus require two degrees of freedom, one of which is the sought for eigenvalue a. Since 
the problem has been linearized, the amplitude of a single asymptotic solution does not 
represent a degree of freedom: one can use the relative amplitudes of two asymptotic 
solutions as the additional parameter to shoot towards the contact line. We thus need to 
identify two linearly independent solutions that satisfy the boundary boundary conditions 
( |4.11|4.12[ ). 



There are two asymptotic solutions of the type: 

zMl ' 



2(a + l)2 ln^{z/c) 
h[ = az"-i ( l + ^ 



Kl 



2(a + l)2 ln\z/c) 
6Ca z"+i 



a + 1 lii%z/c) 



Fi = —— z"+i . (4.13) 
a + I 

These exist for the two roots of 

a'^ + 2a-q^ = a± = ±^/l + q^ - 1 . (4.14) 

Since we want hi to be bounded, only the solution a+ is physically acceptable. In- 
terestingly, the mode hi cx z"+ precisely has the well-known Laplacian structure of 
exp(— + iqy) when transformed in the frame where the bath is horizontal - see Ap- 
pendix [X] This mode thus corresponds to a zero curvature perturbation of a static bath, 
with no liquid flow. Indeed, no flux crosses the bath since Fi ^ 0. 

However, the motion of the contact line implies that liquid is being exchanged with 
the liquid reservoir, so we require an asymptotic solution that has a nonzero value of Fi. 



hi 
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We found that the corresponding mode has the following structure: 

q'Cjz) ^ 1 
z ln^(z/c) 

^, ^ _^ 3_ 

^ zln^(z/c) zln*(z/c) 

ni^q\l + q^)Ciz) 

Fi^^-q^{l + q^) . (4.15) 

C(z) = / dt . (4.16) 

^ ' Jo In^t/c) ^ ' 

This integral can be rewritten in terms of a logarithmic integral using partial integration, 
but this does not yield a simpler expression. 

For completeness, let us also provide the fourth asymptotic solution of this fourth order 
system: 

1 / 3Ca 1 

hi = - [I- 



where C stands for 



z\ l + <?Mn"(z/c) 

' z^\ l + 9Mn2(z/c) 
2^ 3Ca 1 



l + g2 ln'(z/c), 

Fi = crln(z/c) , (4.17) 

which clearly violates the boundary conditions Eq. ( 4.11|4.12 ). 

To summarize, there are two asymptotic solutions that are compatible with the bound- 
ary conditions at the bath. Their relative amplitudes can be adjusted to satisfy one of the 
two boundary conditions at the contact line. The numerical shooting procedure allows 
finding the eigenvalue aq for which also the second boundary condition is obeyed. 



5. The dispersion relation 

5.1. Numerical results 

Let us now discuss the dispersion relation of contact line perturbations obtained within 
the hydrodynamic framework. For fixed microscopic parameters, the relaxation rate 
depends on the capillary number Ca and the dimensionless wave number q that has been 
normalized by the capillary length. This relation will be represented by the function 
CTg (Ca), which has the dimension of the inverse time-scale 7/(77^7). From the definition 
Eq. ( |4.1[ ), it follows that a is positive for stable solutions. 

The dispersion relations are summarized by Fig. |7] displaying for various values of 
Ca. For values well below the critical speed Cac, one finds that the relaxation increases 
with q, in a manner consistent with the quasi-static prediction that a cx \q\ for large q. 
The crossover towards this linear scaling happens around q ~ I, and is thus governed by 
the capillary length. 

Close to the critical point, however, we find two unexpected features. First, it is clear 
from Fig. [7)3 that the linear regime disappears, or lies outside the range of our curves. 
We have not been able to extend the numerical calculation to larger values of q due 
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Figure 7. Dispersion relation obtained numerically, (a) Relaxation rates ct as a function of q, 
in units of 7/(7;/^) and l/l^ respectively. The various curves correspond to values of Ca ranging 
from to 7.5 10~^ increasing by steps of 0.5 10~^, plus Cac. (b) Same but close to the critical cap- 
illary immber: from top to bottom, Ca = 0.0075, Ca = 0.0075687 and Ca ~ Ca^ = 0.00758751. 



to intrinsic instability of the numerical algorithm (the presented curves have arbitrary 
precision). Hence, the crossover value for q, denoted by the inverse wavelength 1/Acut, 
increases dramatically close to the transition. Second, we observe a vanishing relaxation 
rate for the mode q = at Cac, or equivalently a diverging relaxation time. However, 
the rates at finite wavelengths remain nonzero at the transition. This is in contradiction 
with the quasi-static theory, Eq. (2.7), suggesting that aq vanishes at all length scales at 
the transition. 

To characterize these behaviors in more detail, it is convenient to use an empirical 
relation for the numerical curves dOndarquhu 1992"| , 




ctq + cToo 7 . (5.1) 



This form contains the two main features of the dispersion: the relaxation rate for the 
zero mode f7o(Ca), and the prefactor in the linear regime (Too(Ca) = limq_^oo <^q/Q already 



defined in Eq. (2.2 1. The cut-off length Acut then characterizes the cross-over between 



the two regimes. The quasi-static prediction would be that Acut ~ 1 and <7oo oc ctq, see 



Eq. (2.9 1. We have found that Eq. (5.1) provides an excellent fit for all data. Only close 
to the critical point, where the linear regime is no longer observed within our numerical 
range, the values of Acut and CTqo are slightly dependent on the choice for the fit. The 
result for cto is completely independent of this choice. 

Let us first follow the relaxation of the mode q — as a, function of Ca. Figure [8^ 
shows that ctq decreases with Ca, so that the relaxation is effectively slowed down. When 
approaching the critical point, this stable branch actually merges with the first unstable 
branch shown in Fig. ^p, the latter giving negative values for (Tq- As a consequence, 
the relaxation rate has to change sign at Cac, so that ctq = at this point. The graph 
in Fig. shows that this the relaxation time diverges as oc 1/VCac — Ca. As we 
argue below, this behavior is a fingerprint of a saddle-node bifurcation. This scenario is 
repeated when following the higher branches of Fig. |4Jd. Indeed, one finds a succession of 
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Figure 8. (a) Zero mode relaxation rate tro as a function of Ca (solid line). The dashed 
line shows the quasi-static approximation, (b) Same but plotted in Log-Log coordinates as a 
function of Cac — Ca. 



saddle-node bifurcations at which ctq changes sign. In Fig. [8^ this manifests itself as an 
inward spiral, so that the solution with Zc\ — > oo has cr = 0. The dotted curve in Fig. [8^ 



has been obtained from the quasi-static prediction, Eq. (2.5), relating do to the curve 
Zci(Ca) displayed in Fig.jljs: the agreement is excellent. 

This agreement is in striking contrast to the discrepancy at small wavelengths. These 
are represented in Fig. [9^ through (Too(Ca). The comparison with quasi-static theory 
(dotted line), reveals a significant quantitative disagreement for all Ca. However, the 
most striking feature is that CToo diverges near the transition. This suggests that for large 
q the relaxation rates increase faster than linearly, so that the quasi-static theory breaks 
down even qualitatively. A direct consequence is then that Acut — > 0, as can be seen from 
Fig. |9|3. At the critical point Eq. (5.1 1 reduces to dg ~ cr oo Acut 9^/2 for small q, so that 
Acut cx l/cToo- These results underline the quahtative change when approaching Cac. 

5.2. Interpretation 

We propose the following interpretation for the behavior near the critical point. We have 
seen that the q = mode is well described through a standard saddle-node bifurcation, 
which has the normal form 

. ^'''^ 
For positive /i, this equation has two stationary solutions, namely A± — 'i^^/Ji. Linear 

stability analysis around these solutions shows that the A^ solutions are stable while 

the A_ are unstable, and the corresponding relaxation rates scale as cr = ±2,Jfl. So 

indeed, our numerical results for q = Q are described by the saddle-node normal form, 

when taking fx cx Ca — Cac and A — \pl — z^. 

Making an expansion around the critical point that incorporates slow spatial variations 

in y, one would expect the following structure 



dA 



D- 



d^A 



(5.3) 



Due to the symmetry y 



dt dy^ 
-y, the single derivative of y can never emerge. This then 
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Figure 9. (a) Asymptotic relaxation rate aoa = \imq->ao <yq/q as a function of th e cap illary 
number Ca (solid line). The dotted line shows the quasi-steady prediction, Eq. 1 2.1 1, (b) 
Crossover wavelength Acut as a function of capillary number Ca. 



yields a dispersion relation 

cr = <Jo + Dq^ + 0{q'^) . (5.4) 

This explains the observation that for finite q the relaxation rates remain finite at 
the critical point, even though (Tq = 0. When comparing to Eq. (5.11, one finds that 
D = cr ooAcut/2. This value decreases with Ca but remains finite at the transition. Inter- 
estingly, however, the dependence D(Ca) appears to extrapolate to zero only about 1% 
beyond Cac. 



6. Discussion 

We have performed a hydrodynamic calculation of perturbed receding contact lines, in 
which viscous dissipation has been taken into account on all length scales (from molec- 
ular to macroscopic). This goes beyond earlier work by Golestanian & Raphael and 
[Nikolayev fc Beysens 20031 which all dissipation was assumed to be localized at the 
contact line and described by an apparent (macroscopic) contact angle 6a- 

In the first part of the paper we have revealed the bifurcation diagram for straight 
contact lines, which turns out to be much richer than expected from the simplified quasi- 
static approach. Instead of a single saddle-node bifurcation at the critical capillary 
number Cac, we find a discrete series of such bifurcation points converging to a second 
threshold capillary number Ca* (Fig.|4]). Interestingly, the latter solutions have been ob- 
served experimentally as transient states towards liquid deposition QSnoeijer et al. 2006 1 . 
These experiments showed that the wetting transition occurs already at Ca*, and hence 
before the critical value Cac at which stationary menisci cease to exist. Since we have 
found the lower branch of Fig. |4]to be linearly stable at all length scales, this subcritical 
transition has to be mediated by some (unknown) nonlinear mechanism. Let us note 
that similar experiments using thin fibres instead of a plate suggest that it actually is 
possible to approach the critical point (jSedev fc Petrov 199111 . It would be interesting to 
investigate the bifurcation diagram as a function of the fibre radius r, where the present 
work represents the limit r/l^ —f oo. 

The second part concerned the relaxation of perturbed contact lines. At long wave- 
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lenghts, qly ^ 1, we have found that the relaxation obtained in the hydrodynamic 
calculation is very close to the quasi-static prediction. The quasi-static model is based 
upon the equilibrium contact line position at steady-state, as a function of the capillary 
number: it treats the perturbations as a small displacement of the contact line, Azci, that 
induces a change in the contact line velocity ~ dCa/dzd. A positive (negative) derivative 
indicates that the contact line is stable (unstable). This argument does not involve the 
apparent contact angle: it holds as long as the interface profile relaxes adiabatically along 
stationary or steady meniscus solutions. The long wavelength theory therefore relies on a 
"quasi-steady" assumption, and not so much on the interface being nearly at equilibruim 
(quasi-static) . We wish to note that the physics is slightly different for a contact line on a 
horizontal plane, for which there is no equilibrium position due to gravity. For an infinite 
volume translational invariance implies that (Jq = (ISekimoto et al. 1987|l . while drops 
of finite volume has a finite resistance to long wave-length perturbations. This nicely 
illustrates the importance of the outer geometry. 

For small wavelengths, ql^ 1, we found that the quasi-static theory breaks down. 
Away from the critical point we still observe the scaling a oc \q\, as proposed by 
[Joanny fc de Gennes 1984 This scaling reflects the "elasticity" of contact lines, rep- 



resenting axL increase of surface area, and thus of the surface free energy, proportional 
to Quantitatively, however, the quasi-static approximation is not able to capture 
the hydrodynamic calculation. The disagreement becomes even qualitative close to the 
critical point: finite wavelength perturbations do not develop the diverging relaxation 
times predicted by [Golestanian fc Raphael 2001a Also, the scaling a cx \q\ is found to 
cross over to a quadratic scaling a (x q^. 

These results have a clear message: viscous effects have to be treated explicitly when 
describing spatial structures below the capillary length. Namely, the viscous term in 
Eq. (3.3 1 becomes at least comparable to gravity at this scale. Therefore one can 
no longer assume that viscous effects are localized in a narrow zone near the contact 
line: the perturbations become comparable to the size of this viscous regime. However, 
even if contact line variations are slow, a complete description still requires a prediction 
for Zci(Ca), or equivalently 6a{Ga). As was shown by [Eggers 2004[ this relation is not 
geometry-independent so one can never escape the hydrodynamic calculation. 

Our findings provide a detailed experimental test that, on a quantitative level, are 
relatively sensitive to the microscopic physics near at the contact line. In our model we 
have used a simple slip law to release the singularity, but a variety of other mechanisms 
have been proposed previously. The other model parameter is the microscopic contact 
angle 9c\, which we have simply taken constant in our calculations. In a forthcoming paper 
we present experimental results and show to what extent the model is quantitatively 
accurate for the dynamics of contact lines. 
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Appendix A. Perturbation of the static bath away from the plate 

At large distances from the plate, the behavior of the static bath is more conveniently 
described through the function 2surtacc(a^j 2/) — z{x,y). Denoting x positive away from 
the plate, we find asymptotically that z = z' = z" = as a; — > oo. The equation for the 
static interface then simplifies to 

V^z^z, (Al) 
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where we have put 1^ = I. The basic profile is simply exponential 

zo^Ae-^, (A 2) 

while transverse perturbations e*''^ decay along x as 



zi = e-^^^. (A3) 

We can thus write 



i = Ae-^ + ee'«^(e~^)vi+9 , (A 4) 

and compare this to the representation ^surface — h{z,y) 

- ln{z/c) + ee"^y hi{z). (A 5) 



Inserting this x = h{z,y) in Eq. (A4|, and identifying z — z, we obtain to lowest order 
in £ 



Az ( A7\ Vi+s"" 

z= — (l-ee*9''/ii(z))+ee^9M — j , (A 6) 

so that Aj c — 1 and h\(^z) — z", with a — \/T+~q^ — 1. So, the exponential relaxation 
in the frame (a;, y) translates into a power law for hi{z). 

Appendix B. Quasi-static approximation 

In this appendix we derive the quasi-static results summarized in Sec. [2j To perform 
a linear stability analysis, we write the interface profile as 

h{z, y, t) - ho{z) + e h,{z) 6^-*+^??/, (B 1) 

k{z, y, t) = Ko{z) + e (B 2) 

where k is twice the mean curvature of the interface. In the quasi-static approach, the 
basic profile hg^z) can be solved from a balance between capillary forces and gravity, so 
that the scale for interface curvatures is the capillary length = y^j/pg. If we consider 
modulations of the contact line with short wavelengths, 1/q ^ Ij, one can thus locally 
treat the unperturbed profile as a straight wedge, ho^z) — {zd — 2)tan^a, where the 
position of the contact line is denoted by Zc\. Since gravity plays no role at these small 
length scales, one can easily show that the perturbation should have zero curvature, i.e. 
Ki{z) = 0. In the limit of small contact angles, for which we can simply write 

V2(/ii(z)e-"*+*«^) , (B3) 



one directly finds that ki ~ leads to a perturbation decaying exponentially along z, as 
hi{z) — 6^1'"^^='^^^ The length scale of the perturbation is then simply 1/q. This can 



be generalised using the full curvature expression Eq.( 3.4 1: inserting the linearization 



Eq. (B 1 1, and taking dzho = tan 0a, dzzho = 0, one finds 

Kl = (cos 0a)3/2 (dzzhi ^/ii) ■ (B 4) 

V COS^ Oa J 

Hence, the condition that ki = yields: 

/ii(z) = e-l«l^ . (B5) 



From Fig. |3^ it is clear that the "advanced" part of the contact line, with positive Az, 
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has a smaller apparent contact angle than the wedge. The remaining task is to relate the 
quantities s, Az and A6, and to impose the correct boundary condition through ^^(Ca). 
We now introduce the representation 



h{z,y,t) = ho{z - Az) + eh^iz ^ Az) e-''*+"^y , 

in which the position of the contact line is explicitly shifted to z = Zd + Az. 
hi{zci) = 0. Linearizing this equation around z — z^i, this can be written as 



(B6) 
so that 



h{z,y,t) = ho{z) + e \^hiiz)e-^'+"^y - 
= ho{z) + e (hi{z)e-'''+"'y + 



dh{) 


Az\ 


dz 





+ 0{e') 



tan ( 



Az 

e 



+ 0{e^). 



(B7) 



Comparing to Eq. (Bl) with /ii(zci) = 1, one thus finds that to lowest order Az = 
£6^'^*+'^^/ tan Writing dh/dz = — (tan^?^ + Atanfl), one furthermore finds from 
Eq. ( |bT| 

, tan On . , ^ 

\q\^::::^Az. (B8) 



A tan 6* : 



— ut+iqy 



m 

COS Oa COS da 

The final step is to use the empirical relation between 9a and Ca to relate the variation 
in contact angle to a variation in the contact line velocity, Ud ^ U — dAz/dt: 

-1 



dAz 
"df 



1 



ACa: 



.1 



d tan 9n 



V 

\q\j tan.9a f dtanOa 
r] cos 9a \ dCa, 



dCa. 



A tan 6* 



Az 



(B9) 



This indeed results into an exponential relaxation Az oc e with a relaxation rate given 
by Eq. 

For contact angles close to 7r/2, one can analytically solve the crossover from small 
to large wavelengths ( [Nikolayev fc Beysens 2003| . In this case the basic profile is nearly 
flat in the x direction. Characterizing the free surface by -Zsurfaco(a;, 2/), one easily finds 
that perturbations of the contact line decay along x over a distance 1/^^?^ + 1/^^ - see 
Appendix [A} Hence, 

^ ~ ./l + ((7/^)2 for 9a ~Tr/2. (BIO) 



This is consistent with Eq. (2.7l, since g(7r/2) = 1. 
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